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Abstract 

Linked cluster expansions are generalized from an infinite to a finite vol- 
ume. They are performed to 20th order in the expansion parameter to ap- 
proach the critical region from the symmetric phase. A new criterion is pro- 
posed to distinguish 1st from 2nd order transitions within a finite size scaling 
analysis. The criterion applies also to other methods for investigating the 
phase structure such as Monte Carlo simulations. Our computational tools 
are illustrated at the example of scalar O(N) models with four and six-point 
couplings for N = 1 and N = 4 in three dimensions. It is shown how to 
localize the tricritical line in these models. We indicate some further appli- 
cations of our methods to the electroweak transition as well as to models for 
sup er conductivity. 
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1 Introduction 



The phase structure of models for strong and electroweak interactions has been a 
topic of intensive research in the past. In spite of numerous investigations some cen- 
tral questions are still open. To these belong the nature of the chiral/ deconfmement 
transition in QCD for physical values of the current quark masses and the strength 
of the electroweak transition for the physical (so far unknown) Higgs mass. In both 
realms one has to account for nonperturbative coupling regions. Thus it is natural 
to choose the lattice regularized version of these theories to study their phase struc- 
ture. Most applications are performed with Monte Carlo simulations, which are an 
appropriate tool to study the critical region. Monte Carlo simulations are restricted 
to a finite volume. Thorough extrapolations to the infinite volume limit from a finite 
size scaling analysis is in general expensive and sometimes impracticable for lattice 
sizes which are realistic for QCD or for the standard model QXJ] . 

Convergent expansions such as Linked Cluster or Hopping Parameter Expansions 
(HPEs) provide an analytic alternative to Monte Carlo simulations. They may also 
serve as a convenient supplement to numerical calculations. Originally they have 
been developed in the infinite volume. In contrast to generic perturbation theory 
about noninteracting fields, HPEs are convergent power series expansions about 
completely disordered lattice systems. Under certain conditions their convergence 
radius can be directly related to the location of the physical singularity. Hence, 
similarly to Monte Carlo simulations, HPEs can be applied to the phase transition 
(critical) region, if the order in the hopping parameter k is just high enough. Thus 
the transition region is accessible from the high temperature (symmetric) phase. 

Hopping parameter expansions have a long tradition in statistical physics (|2], [3], 
f|] and references therein). Their generalization and application to particle physics 
have been pioneered by Liischer and Weisz . Luscher and Weisz studied a lattice 
$ 4 -theory close to its continuum limit in four dimensions || |7j . Recently, the HPE 
has been generalized to field theories at finite temperature ||. The generalization 
is twofold. First of all one has to implement a toroidal symmetry in one direction of 
finite extension, say L . In this context the temperature T is given by T = Lq 1 in 
lattice units. Second, the highest computed order in the expansion parameter has 
to be increased, because the toroidal (temperature) effect on the critical coupling 
is rather small. Typically, the critical hopping parameter k c changes only by a few 
percent even on a 4 x oo 3 lattice compared to the oo 4 lattice. The graphs of the 
expansion can only "feel the temperature", if they are able to wind around the 
torus in the temperature direction. In general, the largest possible winding number 
should be larger than one to induce a measurable effect. In || the 18th order has 
been used to determine the critical behaviour of the finite temperature $ 4 models 
with O(N) symmetry. Meanwhile, the 20th order of this expansion is available for 
2-point susceptibilitites, see below. 

In this paper we extend the HPE to a finite volume, i.e. to a lattice with toroidal 
symmetry in all directions. We propose criteria to distinguish first from second order 
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transitions (and crossover phenomena) both in an infinite and in a finite volume. It 
is the fate of power series expansions that one cannot work at the singularity k c , one 
can only come close to it, the closer, the higher the order in the expansion. Thus 
we need a criterion that works slightly below k c . As such a criterion we propose 
a so called monotony criterion which is based on the specific volume dependence 
of truncated correlation functions close to but not at the transition point. The 
criterion includes both order parameter susceptibilities and other singular response 
functions such as the specific heat. Decrease or increase with the volume identifies 
first or second order transitions, respectively. Although the monotony criterion has 
been developed in the framework of HPEs, it is not restricted to this case. It can 
be used in other methods as well, in particular in Monte Carlo simulations. 

As a second application of the HPE in a finite volume we calculate an effective 
potential up to 16th order in the hopping parameter. The shape of the effective 
potential further characterizes the type of transition. The coexistence of distinct 
minima at the critical point provides another possibility to calculate k c in a finite 
volume. 

The criteria will be applied to scalar 0(N)-models with $ 4 and $ 6 self- inter- 
actions in three dimensions. These models allow for various first and second order 
transition regions in the bare coupling constant space. For fixed couplings the 
phase transitions will be considered as a function of k. The parameter k may 
be identified with an inverse temperature 1/T of a classical system with the same 
action in three dimensions. Thus we will sometimes replace k c by T~ x . In connection 
with a field theory in four spacetime dimensions the three dimensional model may 
be considered as an effective description of the four dimensional model at finite 
temperature, arising in a process of dimensional reduction. In a four dimensional 
theory at finite temperature one should distinguish between k and T _1 . 

The scalar 0(N)-models contain a number of interesting special cases. If the 
four-and six-point couplings are sent to infinity in an approriate way, we obtain 
"diluted" O(N) models, i.e. Heisenberg models with additional occupation number 
variables, for N = 1 we have a diluted Ising model. The case of N = 4 and 
pure quartic selfinteraction is assumed to share the universality class with QCD in 
the limit of two massless flavours. It also corresponds to the scalar sector of the 
electroweak standard model. A $ 4 + <3> 6 -theory exhibits a tricritical point (line) 
for a fixed (varying) six-point coupling. Such a tricritical point is observed in a 
liquid mixture of He 3 /He A . Recently it has been also proposed as candidate for 
representing the universality class of tricritical QCD |]10| , |i"T| . (Tricritical QCD 
means QCD with vanishing up and down quark masses and a strange quark mass 
which takes a critical value, at which the chiral transition changes its order.) We 
indicate how to localize the tricritical line in a $ 4 + $ 6 -theory with our methods. 

The outline of the paper is as follows. In section 2.1 we summarize the main 
results for HPEs from ||. It basically serves to fix the notation. We then extend 
the HPE in an infinite volume to a graphical expansion in a finite volume (section 
2.2). In section 3 we give two criteria to distinguish 1st from 2nd order transitions: a 
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precise formulation of the monotony criterion (section 3.1), and an effective potential 
evaluated in the HPE in a finite volume (section 3.2). In section 4 we apply these 
criteria to three-dimensional scalar O(N) models with renormalizable interactions. 
To get a first estimate on the phase structure in bare coupling parameter space, we 
study the large coupling limit by a saddle point integration. Another estimate for 
the location at finite couplings is obtained from a hopping-mean-field analysis. This 
approximation amounts to a tree level evaluation of the HPE (Section 4.1). After 
this preliminary study of the phase structure we present a more detailed investigation 
by means of the HPE. In the infinite volume limit, plateaus of critical exponents as 
obtained from the linked cluster series are proposed as criteria to identify the various 
universality classes of the critical region of the theory (Section 4.2). In Section 4.3 
we discuss the finite volume behaviour of various quantities. The shift in volume of 
the critical coupling k c , defined here as the radius of convergence, is compared to 
the scaling behaviour which is expected for the shift of the maximum of the order 
parameter susceptibility. The monotony criterion and the effective potential are 
evaluated for points both in the first and second order coupling region. Finally we 
show how to locate the tricritical line. In section 5 we summarize our results and 
give an outlook to further physical applications. 



2 Hopping parameter expansions for the critical 
region 

2.1 General framework 

Linked cluster expansions provide a convenient tool for both numerical and analytic 
studies of lattice field theories. The typical expansion parameters are the coupling 
strengths between fields at different lattice sites @, [|. In contrast to saddle point 
expansions which are at most asymptotically convergent, series resulting from HPE 



are absolutely convergent for sufficiently small couplings [I2|. In this sense they 
can be viewed as generalized high temperature expansions. If in addition the sign 
of susceptibility series is uniform, the radius of convergence identifies the phase 
transition, i.e. the critical temperature. 

In order to extract quantitative information on the critical behaviour one has to 
get sufficiently close to the critical point. The price to be paid is a computation to 
high orders. The realization of such expansions by convenient algorithms with the 
aid of computers has been pioneered by Luscher and Weisz ||. Recently, progress in 
various ways has been made to extend the length of strong coupling series [§, [13], [l4j . 
Normally, these expansions are set up in an infinite volume. In |J, the techniques 
have been further developed in such a way that the expansions can be reliably applied 
to lattices of non-trivial topology. In particular, it turned out that the highest order 
in the expansion had to be further increased for measuring effects from topology. 
The improved techniques have been applied to scalar O(N) models with quartic 
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interaction on 4- dimensional finite temperature lattices |J. The critical exponents 
could be shown to agree with the critical indices of the corresponding (dimensionally 
reduced) 3-dimensional models. 

In the following we summarize the main formulas from || || to fix the notation 
and to set up the expansion scheme that later will be generalized to a finite volume. 
We consider a D-dimensional hypercubic lattice A = x£ = Q 1 Z/L i , with L ; e N an 
even number or with Lj = oo. Periodic boundary conditions are imposed for each 
finite Li. The restriction to even Li leads to a considerable reduction of the number 
of contributing graphs because it implies that each loop must have an even number 
of lines. The class of models we discuss are described by the partition function 

JV 



Z(J,v) = fl[d N <S>(x) exp(i 

x£A x^tyeA a, 6=1 

■ exp (-J2°S (*(*)) + E E U^a(x)), (1) 



x£A a=l 



where $ denotes a real, iV-component scalar field, J are external sources, v a b(x,y) 

o 

denote the hopping couplings. The ultralocal part of the action S, which depends 
only on one lattice site, is chosen to be O(N) invariant. It should guarantee the 
stability of the partition function (|l|) for sufficiently small v(x,y). Throughout this 
paper we consider as an example the action of a $ 4 + <3> 6 -theory 

= $ 2 + A($ 2 - l) 2 + (j($ 2 - l) 3 , (2) 

which exhibits a phase structure with both first and second order transitions. We 
emphasize, however, that the general techniques are not restricted to this case. 

Fields at different lattice sites interact with the hopping coupling v a t(x,y). For 
the case of nearest neighbour interactions, it reduces to 

, , _ ( 2k 5 a b , x, y nearest neighbour . . 

W,Z/J = | * Q ; otherwise, {6} 

where k is the so called hopping parameter. The nearest neighbour property should 
be understood modulo the torus lengths. Henceforth we consider only nearest neigh- 
bour interactions. 

The generating functional of connected correlation functions is given by 

W(J,v) = \nZ(J,v), 
W^ a Jx u ...,x 2n ) = <$ ai (x 1 )---$a 2 „(a; 2n )> c (4) 

-W{J,v) 



dJ ai (xi) ■ ■ ■ dJ a2n (x 2n ) 



J=0 
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In the following a major role is played by the connected 2-point function and the 
corresponding susceptibility X2 and moments /i 2 , defined according to 

x 

5 a ^2 = < ®a(x)$ b (0) > C . (5) 



x \ i=0 



In field theory, it is convenient to define the renormalized coupling constants via the 
vertex functional 



r(M) = W(J) — ^ J(x) ■ M(o 



xeA 



= E^i E ri 2 i 2n (x 1 ,...,x 2 „)M ai (x 1 )---M a2n (a; 2n ), (6) 

n>0 ai,...,a,2n 

M a {x) = a = l,...,n. 

The standard definitions of the renormalized mass mjj (as inverse correlation length) 
and the wave function renormalization constant Zr are 

f^fo-p) = -^(™ 2 R +p 2 + 0(p 4 ))5 a , b asp^O, (7) 

where ~ denotes the Fourier transform. Eq. ([?]) implies that 

2 

m\ = 2D— , Z R = 2lA (8) 

^2 A*2 

The critical exponents 7, u, rj are defined by the leading singular behaviour at the 
critical point k c , 

lnx2 - -7ln(« c — «), 

lnm^j ~ 2z/ln (k c — k) , as /t k c , (9) 
In ~ vi] In (k c — k) , 

such that urj = 2u — 7. 

If the interaction part (^) of the action is switched off, i.e. v = 0, S($>,v = 0) = 

o o 

S the partition function factorizes, and in turn W(J, v = 0) = W 

(J(x)). In particular, 



W(2«) (~ „ \ _ J f2n-lVI^ 2 «( ai ' • • • ' ° 2n ) ' for X! — X 2 — • • ■ — X 2n 

, otherwise 

(10) 
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with 

d 2n 



v.. 



(11) 

J=0 



2n dJf 

and Cin totally symmetric coefficients in cij, i = 1...2n. 

In practice, the vertex couplings v 2 are obtained from the relation 



n>l 



wi+X^^Gm (12) 



n>l 



with 

v 2n— 5 , UoJ 

/ cF$exp(- S ($)) 

or, alternatively, recursively from Dyson-Schwinger equations. 

The linked cluster expansion for W is the Taylor expansion with respect to v(x, y) 
about this decoupled case, 



W(J,v)= (exp^V^i/) — j A w{J,v) 



(14) 



v=0 



The corresponding expansions of correlation functions are obtained from ( JTjj) by (Q). 
Susceptibilities become power series in k with a nonvanishing radius of convergence. 

The management of such an expansion is conveniently done by means of a graph 
theoretical device. Correlation functions are represented as a sum over equivalence 
classes of graphs, each class being endowed with an appropriate weight. In order to 
make high orders in the expansion feasible it is necessary to introduce more restricted 
subclasses of graphs such as 1-particle irreducible (1PI) graphs, 1-vertex irreducible 
graphs, and renormalized moments. The correlations are then represented in terms 
of the latter two. For further details we refer to || |8| . 

The weight of each graph decomposes into a product of its (inverse) symmetry 
number, the 0(N)-group factor, and the lattice embedding number. It is only the 
latter one that depends on the topology of the particular lattice which is involved. 
In the next section we outline the modifications of the embedding numbers due to 
a finite volume. 



2.2 Extension to the torus 

Let us consider a correlation function, such as in ([5]), on a D-dimensional lattice of 
size L x L% x ... x Lu-i with periodic boundary conditions. Except for a trivial 
volume factor, the embedding number Ir(L , . . . , Ld-i) of a connected graph T 
counts the number of possible ways F can be embedded on the lattice. Embedding 
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means a mapping of every vertex v of T onto a lattice site x(v) = (rr , • • • , xd-i){v) 
consistent with the topology of T. Every two vertices have to be mapped to nearest 
neighbour lattice sites if they are neighboured vertices of T, i.e. if they are connected 
by at least one line. Selflines do not exist. Otherwise the linked cluster expansion 
does not impose any exclusion constraints. In particular, an arbitrary number of 
vertices can occupy the same lattice site. 

It is most convenient to rearrange the computation of embedding numbers in 
terms of random walks 0. Towards this end, the set of vertices of T is devided in 
the disjoint sets of internal 2-vertices and their complement. A vertex v is called 
internal 2-vertex, if it has no external line attached, and there are precisely two 
neighboured vertices of v in T. All internal 2-vertices can be reorganized into so- 
called 2-chains between the remaining vertices in an obvious way. Every 2-chain c has 
an initial vertex i c and a final vertex f c , possibly identical, and it has a length l c > 1, 
where l c — 1 denotes the number of internal 2-vertices of c. Here, for convenience, 
we include l c — 1, in which case c just implies the nearest neighbour constraint on 
i c and f c . On the lattice infinite in all directions the embedding number is then 
written as 

*r(oo*)= e' n^(5^(/ C )(°° c )- ( 15 ) 

(x(v)) c 

The sum runs over all placements of vertices v that are not internal 2-vertices, 
with x(vq) kept fixed for some arbitrary vertex vq to account for the trivial entropy 
factor. The product runs over all 2-chains of T. Af^ y (oo D ) denotes the number of 
free random walks of length I from lattice site x to y. Closed analytic expressions 
can be given for Af^y We notice that 

D-l 

N l £ y {oo D ) ^ only if l-^\ Xi - yi \ > even. (16) 

»=o 

In the finite volume with periodic boundary conditions the topology modifies 



(15) at two places. First, the sites x(v) are now restricted to a cube of size L x 
Li x ... x Ljj-x, and the nearest neighbour constraint, implicit in every 2-chain c of 
length l c = 1 holds modulo the torus lengths. Second, the number of random walks 
M^ y {oo D ) is replaced by 

Ae4(L ,...,^-i) = £ KZ D y+ ,. L (oo D ), (17) 



where 



MOviA'D-lt 



D-l 



L = yjijLj. 



i=0 



The sum in fll7| ) accounts for additional random walks which arise from the possible 
winding around the torus. Due to ([TBI), the sum in ([17]) is finite. 
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3 Finite size scaling analysis with HPE 



In this section we discuss two criteria in the finite volume to determine the order of a 
phase transition. Although the criteria are developed for linked cluster expansions, 
their application is not restricted to series representations. 

The question arises, why one is interested in linked cluster expansion on a torus, 
since the expansion are more easily obtained in the infinite volume. Data of the 
critical region such as the critical temperature are successfully extracted from the 
asymptotic high order behaviour of the coefficients of susceptibility series. The 
typical precision here is within 4-5 digits or even better. 

In general the symmetry of the model alone does not determine the properties 
of a transition. For instance, there may be more than one universality class, cor- 
responding to different ranges in the space of bare actions. As was pointed out in 
0, one should look for plateaus of critical exponents to distinguish between them. 
Problems arise close to the boundary of two such domains. "Smearing effects" oc- 
cur due to the truncation of the series. Both universality domains will influence the 
coefficients, the more the lower the order. This does not pose a problem, as long as 
the domains are sufficiently large and the boundary of the domains extends over a 
negligible coupling range. 

Whereas the location of the phase transition can be determined very precisely 
from the infinite volume series, its order is usually a more intricate question, in 
particular, if the transition is weakly first order. Criteria to distinguish 1st and 2nd 
order transitions can be conveniently worked out in the finite volume. Thus, in our 
case, the finite size effects will be utilized rather than suppressed as artifacts of the 
finite volume. 

A finite size scaling analysis for second order transitions can be based on a 
renormalization group approach, see e.g. [[TJJ. The inverse linear size L' 1 of the 
system is put on an equal footing with other scaling fields like the temperature or an 
external field. The analysis results in a prediction of the leading scaling behaviour 
of a susceptibility x about the critical temperature T c according to 

X (t,L)~|tr P(L/Coo(t)) (18) 

for sufficiently small t and large L, where t is the reduced temperature (T — T c )/T c , 
and 7 is the critical exponent characterizing the divergence at t — 0. The amplitude 
P depends only on L measured in units of the infinite volume correlation length l;^. 
Further properties of P ensure that the height of the peak of the susceptibility at 
T C (L) scales according to x(Pc(L), L) ~ L^^ u \ the width &(L) of the critical region 
shrinks with and T C (L) is shifted compared to T c according to T C (L) — T c 

~ L^ 1 ^ . Here v denotes the critical exponent of the correlation length. 

For a generic first order transition an analogous derivation of the scaling be- 
haviour from first principles is missing in general. The ratio L/^ QO (t) is no longer a 
distinguished scaling variable. In the thermodynamic limit, as T approaches T c , the 
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correlation stays finite and model dependent. The rounding and shifting of thermo- 
dynamic singularities are normally described in a phenomenological approach |16| 



which is based on Monte Carlo results. The height of the peak of the susceptibility 
at T C (L) is expected to scale with L D , both the width a(L) and the shift in T C (L) — T c 
are expected to scale with L~ D as L — > oo, where D denotes the space(time) di- 
mension. A more fundamental finite size scaling theory exists for a class of spin 



models that allow for a particular polymer expansion of the partition function [17 



Whereas the predictions of the finite size scaling of the height of the peak and the 
width of the scaling region reproduce the above mentioned behaviour, the shift of 
the location of the peak is derived to be T C (L) — T c ~ L~ 2D . 

It turns out that our series representations in the hopping parameter k for the 
susceptibilities x cannot be uniquely extrapolated to the critical k c (corresponding 
to T~ as explained in the introduction) to the end that the peak and width of x 
confirm the expected scaling behaviour. However, the specific behaviour of \ close 
to k c is conclusive enough for distinguishing regions of 1st and 2nd order transitions, 
as we will show below, without any need for an extrapolation in k to k c . In addition, 
the scaling of k c (L), defined as the radius of convergence in the finite volume, follows 



the form expected (by analogy) for the shift of the location of the peak from |T7 



The scaling behaviour holds for the Z2 model, but also for the $ 4 + $ 6 -models with 
four components, which are not covered by the analysis of [O], cf. section IV.3. 



3.1 The monotony criterion 

For a certain interval of the scaling region response functions with a nonanalytic 
behaviour in the infinite volume limit show different monotony behaviour for 1st and 
2nd order transitions. Examples for such functions are the specific heat and order 
parameter susceptibilities. They are increasing in volume in a certain neighbourhood 
of T c for 2nd order transitions, and decreasing for 1st order transitions for some range 
in the scaling region, which will be specified below. Since we are not aware of any 
discussion in the literature, although the underlying idea is rather simple, we will 
describe the behaviour in some detail. For definiteness we fix the notation in terms 
of order parameter susceptibilities. 

From the standard finite size scaling analysis one knows that 

x(t,L = 00) < 00 ast-^0 (19) 

for a 1st order transition with a possible discontinuity, whereas 

X (t,L = 00) (20) 

for a 2nd order transition with critical exponent 7 > 0. By definition, regular 
contributions to x ma J be neglected in the scaling region. On the other hand, 
x(0, L) diverges in both crises dbS L approaches infinity. More precisely, at T c , x has 
a "5-function" or power law type of singularity for a 1st or 2nd order transition in 
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the thermodynamic limit, respectively. It is this difference that is responsible for 
the different monotony properties for t 7^ in the finite volume. If t 7^ is small, 
to a given lattice size L s < 00, not too small, one can always find a second size L\ 
with Li > L s such that 

X(t, L>L l )> x(t, L a ) for 2nd order, (21) 
X(t, L>Li)< x(t, L s ) for 1st order, (22) 

that is, x is increasing or decreasing in volume. In the series representation of x 
we can set Li = 00, so that (EHl) and ( P2"|) give strong criteria in the whole scaling 
region where we can use the series. If the volume cannot be made infinite, as in 
Monte Carlo simulations, one is confined to L < L for some L Q . Eq. then still 
holds except for a small neighbourhood of T c , where the susceptibility is increasing 
in volume as in the 2nd order case. The width of this neighbourhood is rapidly 
decreasing with increasing L . 

In the following we make these statements more precise. Let t denote the scaling 
field, i.e. t = (T — T c )/T c , v the inverse volume or an appropriate power of it. The 
infinite volume limit is obtained as v —> from above. The transition range is then 
given by small t and v. Let 

H 2 := {(t,v) G R 2 I v > 0} (23) 

denote a half plane, U C H 2 an open neighbourhood of G H 2 and U* =U\ {0}. 

We discuss the case of a first order transition first. The typical behaviour of a 
susceptibility close to the transition is described by the following 

Definition 3.1 For any uj > we define ^f(JA) as the set of real valued continuous 
functions x '■ ~~ * R- with the following properties. 

1. x e C X (JA* \ {(R, 0)}) ; that is, x is once continuously differentiable for v 7^ 0. 

2. For all nonzero t there is v t > such that for all v < \t\/ut 

I x(t,v) I < u. 

3. With appropriate positive constants c, Ki, K 2 and e we have in U* for v 7^ 



V V 6 

, d , , , K 2 
\j-x{t,v)\ < — . 
ot v z 

As an example, consider the following representation of the magnetic suscepti- 
bility in the volume L D 

X2 (T, L) = cL D exp {-fL 2D {T -T c + dL~ 2D ) 2 ) + r)(T, L), (24) 
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with c, d, e constants and i](-,L) analytic for L < oo, locally uniformly convergent as 
L S oo (so that rj(-, oo) is analytic). With v = L~ D , t — (T — T c )/T c , it is straight 
forward to show that 

X (t,v) := X 2(T,L) 

belongs to ^"(U) for some uo. 

More generally, every function x : W* ^ R of the form 

X(t,v) = I X(t,v) (25) 

belongs to ^ (U) with appropriate u > 0, if the following conditions are satisfied, 
la. x(t,v)eC 1 (W), 

lb. x(t, v) together with its (first) partial derivatives are uniformly bounded in U* . 
2a. / G C 1 (R) is a nonnegative function with /(0) > 0, 
2b. \\m x ^± 00 \x\ 1+e f(x) = for some e > 0, 
2c. (d/dx)f(x) is uniformly bounded on R. 

Any such function has the property that it "approaches 5" locally, i.e. 

lim lim / dt y(t, v) > 
e -0+i;-0+ J_ t AV ' 

and is finite. In this case the limits do not commute. An explicit example for such 
a function is 

f(x) = (— ) x / 2 exp (—cx 2 ) , c > 

71 

with normalization 

/oo 
dx f(x) = 1. 
-oo 

After these preliminaries we now state the following volume behaviour. 

Lemma 3.1 Let w > 0, ^ G (U). There are 5, e > 0, and for every t ^ there 
is u t > such that in U* for all w, v, t with v < 5 and u t w < \t\ < ev, 

X(t,v) > \x(t,w)\. 

In particular the lemma holds for w — 0, i.e. 

X (t,v) > \ X (t,0)\. 

This means that the susceptibility in that part of the transition region where \t\/v < 
e is larger than in the infinite volume limit where w — 0. 
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Proof: Let x G > 0. Differentiability implies that 

X(*> u ) = X(0, u) + t J ds q^x(v, v) 
With appropriate c , K > we have 



r)=st 



X(0,v)>- 
v 



and 



in W*. Hence 



d_ 

di] 



< 



K 



x(t , 0) > si _ m. 



Furthermore, for every t ^ there is f t > such that 

\x(t,w)\ < uj 

for all w < \t\/vf Finally, we choose e = c /(2K) and 5 = c /(4w) and get for v < 5 
and v t w < \t\ < e v 

xM>^- — = f >lx(f,«>)l, 

thus it follows the lemma. □ 

Now we come to the second order transition. In contrast to the first order case, 
at a second order transition the order parameter susceptibility can be divergent in 
the infinite volume limit as the critical temperature is approached. This is described 
by 



Definition 3.2 For any 7 > we denote by ^>1(IA) ^ e se ^ of functions \ '■ 
that are continuous and satisfy the following conditions. 

1. There are constants A, K , e > such that in U* 

I x (t,0) - A\t\~i I < K |t|- 7+e . 
Furthermore, with appropriate v , C > 0, we have whenever \t\ > v v, 

X (t,v) > C X (t,0). 

2. There are constants rj, B > such that for \t\ < rj v 

\ X (t,v) I < Bv~\ 



R 
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The specific property for a 2nd order transition that the singular part of the 
free energy density behaves as a generalized homogeneous function implies for the 
susceptibility in a volume L D a typical form like 

X 2(T, L) = \T- T c p Q((T - T C )L^) + V (T, L), (26) 

with some 7 > 0. Here 1](-,L) has similar analyticity properties as in ( f24"D above, 
v > is the critical exponent of the correlation length 

Z~\T-T C \- V , (27) 

Q is continuous and behaves as 

lim |x|" 7 Q(a;) = K > 0, 

lim Q(x) = C>0. (28) 

x— >±oo 

The first equation expresses the absence of a nonanalyticity of X2 for finite L, the 
second one its presence in the infinite volume case. With t = (T — T c )/T c , v = L~ x l v 
and 

X(M) := X2(T,L) 

we see that x belongs to ^(^O fo r some 7- 

More generally, every function x : W* — > R of the form 

X(*,«) = \f£) + x(t,v) (29) 
with 7 > belongs to some ^{U), if the following conditions are satisfied. 
1. x(t,v)eC 1 (W). 

2a. / G C\B), and /(0) > 0, 

2b. lim^^-i-oo |x| 7 /(x) = C for some finite C > 0. 



Compared to the first order case (25) the essential difference comes from property 
(2b). 

As an example, 

X (t,v) = {t 2 + v 2 )- {m/2) ,m>0, 

belongs to the class 

For these functions, we have in contrast to Lemma ( pj.l| ) 

Lemma 3.2 Let 7 > and x G ^(U). There are constants v, e > 0, v < e, such 
that for all t, v, w with vw < \t\ < ev 
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The inequality is always true if w = 0, i.e. the susceptibility is always smaller 
than in the infinite volume limit as long as we are in the critical region \t\/v < e. 

Proof: Let \ ^2(^)5 7 > 0. There are numbers C,V, v > such that in IA* 
for \t\ > v w 

x(t,w) > e x (t,o) > v \t\-\ 

Furthermore, there are rj,B > such that for \t\ < 77 v 

\ X (t,v)\ < Ev-\ 

We choose e = min (rj, (V/B) 1 ^ 1 ) and get for v w < \t\ < e v 

( e \ 7 Be 1 
\x{t,v)\ < Bv 7 < B[ — ) < —x(t,w) < x(t,w). 



t / V 



This proves the lemma. □ 



Figure 1: (t,f)-plane for susceptibilities x{^i v ) > in the vicinity of a phase 
transition at (t = 0, v = 0), t denotes the scaling field t — (T — T c )/T c , v is inverse 
to some power of the volume L x with some x > 0. For a 1st order transition, 
x(t,v\) > x(t,f2), whereas for a 2nd order transition x(t, v i) < x{^i v 2)- For the 
shaded part no prediction is made. 




t 



3.2 The effective potential 

Let us briefly discuss another method to determine the nature of a transition that 
will be of use later on. For definiteness we come back to the iV-component scalar 
model as described in Section p]. A possible way to define an effective potential is 
by 

V-V eS (M) = -r(M)| A/=const , (30) 

where V denotes the volume, and T(M) is defined by @. In the symmetric phase 
and in the infinite volume limit, V e s has to be convex. In practice, the right hand side 



of (|30D is obtained as an expansion about M = 0. In the linked cluster expansion the 
coefficients can be expressed in terms of 1PI susceptibilities Xn PI - They are obtained 
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as series representation in k of the truncated susceptibilities by keeping only those 
graphs that are IPI [H. Up to a constant we obtain 



Kff(M) = I 1 - ^ V - 1 -^(M 2 ) 2 (31) 



2 x ipi 4 , (x ipi,4 



2 

1PIN2 



' 1 ^ PI - 1 %r^)(^ 2 ) 3 + 0(M«). 



6! (xD 6 V A6 xf l 

Any nonconvex shape of V e g in the infinite volume limit must be an artifact of the 
approximation scheme. In a finite volume, however, a nonconvex shape in the sym- 
metric phase signals a first order transition, whereas a convex shape is compatible 
with a second order transition. An estimate of the critical coupling is then obtained 
by a root of the coefficient of M 2 in the 2nd order case, and by degenerate values of 
V e g at the trivial and nontrivial minima in the 1st order case. In our applications 
(see Section 4) we have calculated the lPI-susceptibilities xln 1 U P t° 16th order in 
k in a finite volume for n = 1,2, 3. 

4 Applications to scalar O(N) models with <3> 4 and 
<l> 6 -point couplings 

In this section we apply the methods discussed in the previous sections to the three- 
dimensional O(N) symmetric scalar model with N = 1 and N = 4. The model is 
described on a lattice A by the partition function 

Z(k, A, a) = IH d N $(x) exp (2k £ $(x) • ${y) - £ °S ($(x), A, a)), (32) 

xeA x,y NN x 

where the first sum of the exponential runs over unordered pairs of nearest neighbour 

o 

lattice sites, and the ultralocal part S is given by 

S(<S>,\,a) = $ 2 + A($ 2 - l) 2 + a($ 2 - l) 3 (33) 

with a > or a = and A > 0. In contrast to the pure quartic interaction which 
only admits 2nd order transitions, the action (|33| ) allows a richer phase structure 



with regions of 1st and 2nd order transitions due to the additional $ interaction. 
The case of A = 3cr corresponds to a pure $ 6 -theory, whereas A < 3a implies a 
negative quartic coupling. 

4.1 Preliminaries 

To get a first estimate of the phase structure we consider the case of large couplings 
A and a. For finite coupling constants we invoke a hopping mean field analysis. 
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The large coupling limit. To study the limit of large couplings, we set A = aa 
and send o to infinity with a and k kept finite and fixed. This limit is discussed 
for more general contact terms in Appendix |A]. The discussion is based on a saddle 
point integration. As a result we obtain the following behaviour in dependence on a. 
For a > 1 and a < —3 we obtain O(N) Heisenberg models (Ising model for N — 1). 
The range of —3 < a < 1 leads to complete disordering with no phase transition at 
any finite k. The cases of a = —1 and a = 3 are peculiar. The resulting actions 
describe "diluted" O(N) models, with particular values of the couplings. If the large 
coupling limit is performed term by term in the HPE series it can be shown that at 
least for N > 2 the resulting actions again belong to the unversality class of O(N) 
Heisenberg models. 

A hopping-mean-field analysis. To get a first estimate of the phase structure 
at finite couplings, it is instructive to start with a mean field analysis. Together 
with the convexity of the exponential and the positivity of the measure this ansatz 
leads to a complete factorization of the partition function. The hopping mean field 

O 

estimate for the free energy is then derived as follows. Let us define x by 
expx(H) = I d N <$> exp(-^ S ($,A,<r) +H- $), 

J X 

and the expectation value < F >h of an observable F according to 
<F($)> H = exp(-\A\°x(H)) 



I \\ f d N $(x) exp (- J2 S ($) + H ■ $) j 



Here |A| denotes the lattice volume, and H e M. N is an auxiliary field. Note that x 

o 

defined in this way agrees with W as introduced in Sect. y. In particular, for every 
integer n, 

° V 2n (A,tT) = X{H) 

Un l H=Q 

(cp. to (0)). For simplicity, where no confusion can arise, we only indicate the 
dependence on H. We get 



with 



Z(k,\,<t) = < exp (2k Hx)-$(y))> H 

x,y NN 

> exp<2/t ^2 ®(%)-®(y)>H 

x,y NN 

= ex P (-|A|7(#)) 
f(H) = - (x(H) + 6k(V h x(H)) 2 - H -V H £(#)). 
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An upper bound on the true free energy density / is thus given by 

/ < Mf{H). (34) 
A vanishing Hq = is always a solution of the corresponding mean field equation 

d H J(H ) = -d 2 H °x(H ) (l2Kd H x(H ) - tf ) = 0, 
where the derivative is in the direction of H, and H = is a local minimum of /, if 

d 2 H f = d 2 H °x (H ) (l - 12k d 2 H °x (0)) > 0. (35) 

In particular, for small hopping parameter k the model is always in the symmetric 
phase. 

This form of mean field analysis is identical to the tree level of the hopping 
parameter expansion. It can be easily shown that the Lebowitz inequality 

d? H x (H) < for all H > 

together with d 2 H x (H) > for all H ensures that H = is the absolute minimum 
of /. The Lebowitz inequality holds in any case for A > 3cr. Equality in (RS|) along 
with 

47(0) > (36) 
locates a 2nd order phase transition to the spontaneoulsy broken phase at 

K c (A,cr) = 55- . 

12^ 2 (A, a) 

Tricritical points are identified by an equality in (|36"D and 

^7(0) > 0, (37) 

i.e. 

O C o C 

v 4 (A, a) = and f 6 (A, a) < 0. 

We notice that these conditions imply a vanishing two and four-point coupling 
in the effective potential V e ff, Eq.(|3"TD, evaluated to tree level in the HPE, i.e. in the 
hopping mean field approximation. Thus the criterion for tricriticality reduces to 
the familiar one. If we had used the classical potential (^) instead, this would lead 
to a location of the tricritical line in the bare coupling constant space at A = 3a. 

In Table 1 we have listed some results for the 0(4)-model on the location of the 
tricritical line in the (A, cr)-space for several values of a. 

The tricritical exponents in hopping-mean-field establish the results of the most 
naive mean field analysis with the effective potential replaced by the classical po- 
tential. They are a = 1/2,(3 = 1/4,7 = 1,S = 5,u = 1/2 . From the Ginzburg 
criterion one may expect that the only chance where a mean field type of analysis 
may lead to reliable predictions of the singular behaviour in three dimensions is at 
tricriticality. In fact, the susceptibility comes out as volume independent along the 
tricritical line (cf. Section 4.3 below) when it is determined by the HPE analysis. A 
mean field analysis is volume independent by construction. 
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Table 1: Mean field tricritical line for the 0(4) model on the three-dimensional 

qC qC qC 

hypercubic lattice. The line is defined by v 4 — and v e < 0, cx is defined by a =v 6 
/(5(C) 3 )- 



a 


A = aa 


a 


a 


1/1.05 


128.199 


134.609 


-0.7245305 


1/1.1 


53.040 


58.344 


-0.6875585 


1/1.2 


21.822 


26.286 


-0.5920738 


1/1.3 


13.252 


17.228 


-0.5193738 


1/1.5 


7.315 


10.973 


-0.4451270 


1/2.0 


3.401 


6.802 


-0.3839411 


1/5.0 


0.788 


3.941 


-0.3407590 


1/10.0 


0.344 


3.494 


-0.3345469 


1/100.0 


0.031 


3.081 


-0.3306763 


-1/10.0 


-0.273 


2.731 


-0.3284753 


-1/1.1 


-1.391 


1.530 


-0.3519023 


-1.0 


-1.472 


1.472 


-0.3581768 


-2.0 


-2.948 


1.474 


-0.5106579 


-3/1.2 


-7.561 


3.024 


-0.6827047 


-3/1.1 


-17.678 


6.482 


-0.7259408 



4.2 Infinite volume analysis at finite couplings 

So far we have studied the phase structure in the large coupling limit and in a mean 
field analysis for finite couplings. Next we utilize the linked cluster expansions for 
a more thorough study. Susceptibilities are represented as convergent power series 
in the hopping parameter, such as the 2n-point functions 

X2n(K,\,a) = J2tf n) (^°)^i (38) 

and similarly for weighted correlations. These series have been computed to 20th 
order in k for n = 1, to 18th order for n = 2, and to 16th order for n = 3, both in 
a finite and infinite volume. In the infinite volume, the coefficients of the series we 
have explicitly calculated are of equal sign for each series. Under the assumption 
that this behaviour continues to all orders in k, the radius of convergence n c (\, a) 
is identified with the singularity closest to the origin on the positive real axis, hence 
with the physical singularity at the phase transition, independently of the order of 
the transition. 

Well developed methods are known to obtain critical data from the high order 
coefficients of high temperature series J|, |i~9fl . The critical point k c (A, a) is identified 
by the ratio criterion, applied to the coefficients a^ n \\, a). The best choice is the 
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2-point susceptibility, because its series is available to the highest order. The obliged 
regression towards large \i is done according to 



7 V 



,(2) 



(2) 
%-l 



' f l + ^- + 0(^)\ (39) 



with u 2 > oji > 0, ci as fit parameter chosen according to the best x 2 /df fit. 
This procedure is eventually supplemented by a shift of the weak antiferromagnetic 
singularity at — k c to — oo, an improved estimator fit and other known techniques like 
Pade' methods. We know that u)\ = 1 for a leading pole or branch point singularity 
on the real axis, i.e. for 

X2 — A (k c — k)~ 7 as K y K c 

with 7 > and 7^1, and u\ > 1 for 7 = 1. For a 2nd order transition, an 
alternative way to determine the critical point is given by the smallest real solution 
of 

12k cX I p \k c ) = 1, 

as proposed in |§. This condition is equivalent to the identification of the phase 
transition as a zero of the quadratic coefficient of the effective potential, cf. (|31|). 
It turns out to be the most convenient way to determine the radius of convergence 
leading to the highest precision in k c . The lowest precision obtained in this way lies 
within 4-5 digits. Once k c is determined we obtain the critical exponent 7 from 

1 + v(K c r„ - l) = 7 + + (/i-n (40) 

/i 

with 003 and C2 as fit parameters. In a similar way, the critical exponent v is obtained 
by replacing the series of \2 by that of m^ 2 , cf. (||). 

A measurement of critical exponents like 7, u, r] leads to a qualitative plot of the 
phase structure in the(A, <j)-halfplane as shown in Fig. ^. The solid line represents 
the boundary A = A((cr) between the 2nd and 1st order region. To the left of 
this tricritical line the phase transition is of 2nd order. Here, except for the origin 
A = a = 0, we obtain one universality class for every iV with plateaus of critical 
exponents, with values of the iV-component Heisenberg model (cf. e.g. for a recent 
list of those exponents). It is remarkable that this range considerably extends the 
"Lebowitz domain" A > 3cr, where the action is convex and therefore 7 is not less 
than 1 20]. In particular, it includes the full range of A < —3a. In passing we 



mention that for A > the presence of a small nonvanishing $ interaction, that 
is a > 0, considerably accelerates the convergence of the high temperature series 
compared to the case a = 0, i.e. the non-universal remainder of ( |39|) becomes 
smaller. 

Except for the values of the critical exponents, the phase structure qualitatvely 
confirms the mean field analysis of the last subsection. As the tricritical line A = 
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Figure 2: Qualitative plot of the phase structure of 0(N) lattice models in 3 
dimensions. The dashed curves give the lines A = aa with a = 1 and a = —3. The 
solid curve represents the tri critical line At(er). To the left of it the phase transition 
is of 2nd order, to the right of it 1st order. 



A 




a = — 3 



At(cr) is approached, the exponents 7 and v drop continuously from their values in 
the Heisenberg models to the Gaussian values, where the tricritical line is crossed, 
and further to zero. The smooth interpolation between the different exponents is an 
artifact of the truncation of the power series expansions at high, but finite order in 
Kj. It was already pointed out in || that various universality classes lead to smearing 
effects at finite order due to an "interference" of various universality domains. The 
most pronounced plateau structure is obtained for urj, the critical exponent of the 
wave function renormalization constant Zr 

In Fig. |3| we show the results on the exponent vr) for the 0(4) model, obtained 
from the 20th order susceptibility series of X2 and fi2, @, for various a along the 
ray A = (1/2)<t. There is a well established plateau at the left part of the plot corre- 
sponding to the universality class of the Heisenberg model. The plateau at the right 
part is compatible with a range of 1st order transitions, all exponents 7, z/, 77 vanish 
within the error bars. Hence they are compatible with a finite correlation length 
at k c . The stability of the extrapolated convergence radius k c under a variation of 
the truncation of the series suggests that there is really a 1st order transition rather 
than a mere crossover phenomenon. We would like to identify the left boundary of 
this plateau as the tricritical point. This gives us an estimate of about <j t — 9.0. The 
indicated errors in Fig. ^| are obtained as discussed in connection with (fK]). The 
smearing effect does not allow for a more precise location of the tricritical point. To 
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Figure 3: Critical exponent vr\ of the wave function renormalization constant Zr 
for the 3d 0(4) model, obtained from the 20th order susceptibility series for various 
a along the ray A = (l/2)o\ The tricritical point is at a t ~ 9.0. The mean field 
estimate gives a t = 6.8. 
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get a clearer identification of the 1st order transition region and a better localization 
of the tricritical point it is natural to perform a finite size scaling analysis, which is 
the topic of the next section. 



4.3 Results of the finite size scaling analysis 

In Sect. |] we have formulated monotony criteria for response functions \ in the 
scaling region. In absolute value, increase in volume implies a 2nd order transition, 
decrease in volume a 1st order one. The response functions have to be calculated 
at some k close to but not at the critical point k c for two volumes. The value of k 
should be chosen sufficiently close to the transition point to satisfy the conditions of 
the monotony criterion, and sufficiently apart from k c to allow a use of the truncated 
series representation. A choice of k = 0.98k c (A, a) fulfills both restrictions. 

The volumes should be sufficiently large in lattice units to guarantee the appli- 
cability of the finite size scaling ansatz for x- Beyond this generic condition the 
following restrictions arise from the monotony criterion. The smaller one of the two 
volumes should satisfy L x \k — k c \ < 1 with x = 1/u or x = 3, which implies an 
upper bound on L for given k. In practice we have chosen this L between 4 and 12. 
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In the context of the HPE the larger one of the two volumes may be set to infinity. 
The advantage of this choice is that in the 1st order case, decrease in volume holds 
all over the scaling region. In Monte Carlo simulations the second volume is neces- 
sarily finite. In this case L should be large enough so that k lies outside the small 
neighbourhood of k c where x 1S increasing in the 1st order well. The critical 

point K C (X, a), which enters the inequalities on L and k, is determined as the radius 
of convergence in the infinite volume as described above. 

We use the 2-point susceptibility , because its series is available to the highest 
order in k. For given couplings A and a 



M 



xi M) (k, A, a- L) = af (A, a; L) W (41) 

At=0 

denotes the 2-point susceptibility truncated at order M. With 

r M {X,a;L) := 1 - -, (42) 

we know that r^-; L) > for 2nd order and r OC) (-; L) < for 1st order transitions, 
if L lies within the bounds as explained above. The convergence of the series ( pTT| ) 
as M — > oo ensures the same behaviour for finite, but sufficiently large M. 

Fig. H shows an application of the monotony criteria to the 1-component model. 
We have plotted the volume dependence of the ratios 7"m(A, o",L) for a truncation 
at order M = 20 and for various lattice sizes L = (L1L2L3) 1 / 3 , at two points of 
the bare coupling constant space. One point is well inside the 1st order region, the 
other one well inside the 2nd order part of the phase diagram, the different areas 
have been identified by the infinite volume series as discussed in the last subsection. 
Clearly the sign of r2o(-', L) is different in both regions of phase space. It is positive 
for the 2nd order transition and negative for the 1st order transition. The approach 
to the infinite volume limit where rjvf (•; 00) = is fast. 

Next we want to demonstrate how one can utilize the finite volume criteria with 
HPE to get a better localization of the tricritical region. Let us consider the 0(4) 
model and determine the behaviour of the ratios ru along the line A = (1/2)<t. 
From the infinite volume analysis of the last subsection we had obtained o t — 9.0 
as an estimate for the tricritical coupling a t . 

Since r^(A, a; L) > for 2nd order and < for 1st order transitions, the tricrit- 
ical point should be localized at the zero of r M (\,a;L) between these two ranges 
(with nonvanishing slope, i.e. (9rAf (a/2, a; L)/da 7^ 0), suitably extrapolated to 
M,L -> 00. 

Fig. § shows the ratios r M (a/2, cr; L) as function of o for L = 4 and various M 
between (corresponding to the mean field approximation) and 20. The intersec- 
tions of the curves with the r = 0-axis lie in the range of 8 < a < 8.5. The zeroes 
ctm(L), defined by rM(c r M(L)/2,aM(L); L) = 0, depend on the order M at which 
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Figure 4: Dependence of the ratios r 2 o(X,a;L) as defined by fl4"2"| ) on L = 
(LiL 2 L 3 Y^ 3 , at the example of the Z 2 -model. Two points of the phase space have 
been chosen. 

1. (+) A = 15/1.1 and a = 15. The transition is 1st order, r 2 o < 0. 

2. (□) A = 3 and a — 1. The transition is 2nd order, r 2 o > 0. 
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the susceptibility series of \i have been truncated and on the lattice size L. The 
dependence on L for fixed M of tm(A, <j; L) is shown in Fig. || 

Thus a final localization of the tricritical coupling at needs an extrapolation in 
L and M to infinity. Clearly both extrapolations are not independent of each other. 
We should expect that a comparision of two ratios is sensible, i.e. 

r M (-,L)~r M ,(-,L') (43) 

if lattice sizes L, L' and truncation M, M' satisfy 

M _M' 
T ~ U' 

The reason is that M/L is the maximal number of times a graph contributing to 
the series of X2 ("! L) can wind around the volume. Eq. fl44]) then ensures that the 
remaining L-dependence becomes independent of M for sufficiently large M. 

Fig. ^ shows the data (Jm{L) obtained on a 4 3 -lattice (L — 4), and for 11 < M < 
20, as a function of 1/M. The curves show the regression 

6(M min ) 



(44) 



(T M (4) = a t (M min} L = 4) + 



M 
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Figure 5: The ratios r^A, a\ L) as function of a for A = a/2, and L — 4, i.e. on a 
4 x 4 x 4-lattice. 
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as obtained for M min < M < 20, and for M min = 14, 16, 18. In Fig. |S] we show the 
second regression for the resulting a t (M min L / '4, L) in 1/M min , for L = 4, 6 

L e 
er t (M TOi „— , L) = cr t + , 

leading to the final prediction of a t - We have scaled the adjusted M min according 
to (PI). Note that both data sets, on the 4 3 and on the 6 3 -lattice, fall on the same 
straight line within the error bars. We obtain 

a t = 9.454(49). (45) 

In passing we remark that the accuracy of at has been increased by at least one 
order of magnitude compared to the infinite volume analysis (Fig. |). It should be 
noticed that the volume independence of \2 in the scaling region of the tricritical 
point just confirms the validity of a mean field analysis of tricritical exponents. A 
mean field analysis is volume independent by construction. The expected volume 
independence of \2 is confirmed within the error bars. 

In the remainder of this section we discuss the volume dependence of the radius 
of convergence and of the effective potential. 

Shift and scaling of n c (X,a; L). We define K c (X,a;L) in the finite and infinite 
volume as the radius of convergence of the susceptibility series, in particular of X2, 
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Figure 6: Dependence of the ratios r 19 (cr/2, a\ L) as defined by ( |4"2"D on a for various 
volumes. 
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which is known including the 20th order. Some data for k c (L) are listed in Table 
for a three-dimensional 0(1) and 0(4) model. The couplings (A, a) for each N 
(N =1,4) have been chosen deeply in the 1st order transition region. 

From a finite size scaling analysis one expects (cf. Sect. |||) that the data should 
fit with a regression in L according to 

\k c (-; L) - k c (-] oo)| ~ lnc - y T hiL 

for large L, with some constant c and a critical exponent jjt- For the Z 2 model we 
obtain in the 1st order transition region (A = 15/1.1, a — 15.0) according to Table |2| 



In c 

Vt 



4.57(57) 
6.21(32) 



with x 2 /df =0.025. Thus the scaling behaviour is consistent with yx = 2D = 6. It 
confirms the behaviour which has been predicted for the shift of the critical coupling 
determined as the maximum of the susceptibility in a class of models which cover 
the Z2 model |]I5| . Note that it is in disagreement with the Gaussian two-peak model 
[ 16| , predicting a leading finite size correction proportional to L~ D , which one might 
have expected. 

The same scaling behaviour is found for the 0(4) model in the 1st order transition 
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Figure 7: The solution <tm(L) of the equation r M (a M (L)/2,<jM(.L)] L) = for 
L — 4, plotted against 1/M, M is the order of truncation of the suceptibility series. 
Regression is shown as a linear function of 1/M for M > M m i n with M m i n = 
14,16,18. 
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region. Here we get (Table |J) 

lnc = 4.57(120) 
y T = 5.55(59). 

Thus a leading correction proportional to L~ 3 lies clearly outside the error bars. 
This result is remarkable, because the 0(4)-model (which is a Heisenberg model in 
the large coupling limit) is not covered by the analysis of Borgs, Imbrie, Kotecky 
and Miracle-Sole |T7| , |T8|j . Hence the vanishing of the coefficient of the linear term 
in a large volume expansion of the susceptibilities in powers of 1/L 3 seems to be a 
universal feature of a large class of models. 

A measurement of the scaling behaviour of k c (L) in the 2nd order region of the 
0(4) model was not conclusive, because the shift of k c occurs in the 4th or 5th 
digits, hence the finite size effect is hidden in the error. In the 2nd order region of 
the Z2 model we measure a scaling, which is best fitted by an ansatz (cf. Table ^) 

k c (-;L) = « c (-;oo) + — 

with y T = 2/u, and 

lnc = -2.25(129) 
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Figure 8: The "scaled" solutions a t (M min L/4:, L), plotted against 1/M min , for L = 4 
and L = 6. Within the error bars, the data are on a straight line, thus confirming 
the assumptions made in connection with Eqns. (53) and (f4*l|). Linear regression in 
1/M min gives the tricritical point. 




0.12 



y T = 2.83(59) 

with X2/df =0.1. From a renormalization group analysis one expects a leading scaling 
correction proportional to L~ x l v ' . Also here the finite size effect is not large compared 
to the error in k c , thus the disagreement with an L -1//l/ -behaviour may be due to 
the errors in k c and the uncertainty in the involved extrapolations. 

We conclude the discussion of K C (L) with a conjecture concerning the monotony 
behaviour. 

Monotony in k c (L). The results of Table |2| for K C (L) exhibit a further character- 
istic distinction between 1st and 2nd order transitions. Since we are not aware of a 
general proof, we leave it 

Conjecture: The convergence radius k c (L) of the series expansions is monoton- 
ically decreasing with L for 2nd order transitions, and monotonically increasing for 
1st order transitions. 



The effective potential as function of L. Similarly we have measured V e ff, 
Eq. (|3~1~D, for two points well inside the supposed 1st and 2nd order transition regions 
of the 0(4) model (A = 6.0, a = 12.0, and A = 0.90, a = 1.0, respectively). The 
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Table 2: Radius of convergence k c (L) of the HPE series for the 1st order region. It 
is determined from the 2-point suceptibility series, as described in Section 4.2. Data 
are given for the 3d Z 2 and 0(4) model, for various volumes L 3 . 



L 


Z 2 


0(4) 


A = 13.64 , a = 15.0 


A = 6.0 , = 12.0 


K C (L) 


K C (L) 


00 


0.51047(1) 


0.84462(10) 


4 


0.49381 


0.79531 


6 


0.50886 


0.83851 


8 


0.51025 


0.84416 


10 


0.51047 


0.84436 


12 


0.51047 


0.84447 



Table 3: Radius of convergence k c (L) of the HPE series for the 2nd order region in 
the Z 2 model for various volumes L 3 . 





A = 3.0 , a = 1.0 


L 


K C (L) 


00 


0.17316(4) 


6 


0.17393 


8 


0.17339 


10 


0.17327 


12 


0.17327 



xln 1 entering Eq. (EH]) have been evaluated to 16th order in k for n — 1,2,3 in 
a finite and infinite volume. For the 2nd order point A = 0.90, o = 1.0, and for 
volumes L 3 with L = 4, 6 and infinity, V e ff is convex in the symmetric phase up to 
a resolution of 10~ 8 . 

In the first order case (Fig. |^) V e ff is nonconvex in the symmetric phase with a 
barrier height decreasing with increasing L. The three curves correspond to a 4 3 , 6 3 
and oo 3 lattice, respectively. Note that the nonconvex shape which is even seen for 
the oo 3 lattice (lowest barrier in Fig. |S|) must be attributed to the approximation 
scheme, i.e. to the truncation at order M = 16. (As mentioned above, the truncation 
of the series expansion acts similarly to a finite volume cut-off. The finite volume 
leads to a nonconvex shape only in the 1st order case. Thus it is not surprising that 
we find a convex shape in the symmetric phase for all volumes in the 2nd order case, 
in spite of the same truncation in M.) 
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Figure 9: Volume dependence of the effective potential Kff($) on lattices with 
varying size L 3 . The barrier height decreases with L. The curves are obtained for 
L = 4, 6, oo. The parameters are A = 6.0, a = 12.0 in the 3d 0(4)-model. 
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The coexistence of minima leading to the same value of V e ff defines a critical 
coupling k c (L) which need not agree with the finite volume convergence radius k c (L), 
unless the truncated expansions of V^($) in $ and of xln 1 * n K are extrapolated 
to infinity. For parameters chosen as in Fig. |9| we find k c (4) = 0.8311, k c (6) = 
0.8475, k c (oo) = 0.8488, in contrast to « c (4) = 0.79531, k c (6) = 0.83851, k c (oo) = 
0.84462, cf. Table |2], determined as the radius of convergence of the series expansions, 
extrapolated to infinite M. 

Tricritical parameters from V e ff. Vanishing coefficients of the <£> 2 and $ 4 -terms 
associated with a qualitative change in the shape of V e ff provide a further possibility 
for localizing the tricritical couplings. The analytical dependence of the coefficients 
on A and o is rather indirect. It is easier to determine At, Ot by the first occur- 
rence of a nonconvex shape of V e ff in the symmetric phase, coming from the 2nd 
order transition region. The highest order, which is so far available for the six-point 
susceptibility \\ P1 \ entering the $ 6 -coefficient of V e ff, is 16. The number of con- 
tributing graphs to order 16 is comparable to the number of graphs contributing to 
the two-point susceptibility to order 20 and the four-point susceptibility to order 
18. An inclusion of higher powers in $, say <3> 8 , <3> 10 -terms, would further reduce the 
order in k which is tractable. Thus we retain from further extrapolations, but give 
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bounds on a t , derived from V e //($) to 0(/t 16 ). They are 

9.75 <a t < 10.0. (46) 

The resolution in V e ff, within which no nonconvex shape was seen up to a = 9.75, 
was 10~ 7 . Compared to the more reliable result of Eq. (f45"|), the evaluation of V e ff 
seems to lead to an upper bound on cr t , given by Eq. (PH). 



5 Summary and Outlook 

In this paper we have generalized hopping parameter expansions from an infinite to a 
finite volume. The combination of performing the expansions in a finite volume and 
to a high (20th) order in the expansion parameter has turned out as a useful compu- 
tational technique to approach the critical region from the symmetric phase and, in 
addition, to characterize the type of transition. First and second order transitions 
have been distinguished by various criteria: the monotony criterion refering to the 
L-dependence of response functions (here illustrated at order parameter suscepti- 
bilities), the scaling and monotony of the radius of convergence k c (L) as function of 
the linear lattice size L, and the effective potential as function of L. In particular, 
it is the different L-dependence of the order parameter susceptibility for 1st and 
2nd order transitions in the scaling region (close to but not at the transition point) 
which allows us to localize tricritical points. The monotony criterion can be applied 
to Monte Carlo simulations as well, the involved two volumes should be sufficiently 
large, but both may be finite. 

We have applied these methods to renormalizable O(N) models in three dimen- 
sions. The plateau structure in the critical exponent vr\ in the infinite volume has 
revealed two universality classes belonging to an 0(iV)-Heisenberg model and to a 
Gaussian model. Apart from the "trivial" Gaussian behaviour at vanishing four and 
sixpoint couplings, we get Gaussian exponents along a tricritical line separating 1st 
and 2nd order domains. The existence of the 1st order domain and the tricritical 
line is based on the presence of the <3> 6 self- interaction. The O(N) symmetry alone 
does not determine the critical behaviour. 

Several extensions are at hand. The first one is from 3 dimensions to field 
theories in four dimensions at finite temperature. For a $ 4 + $ 6 -type theory in 
four dimensions at finite temperature we expect qualitatively the same infrared 
behaviour and phase structure as for three dimensions, but different values for the 
critical couplings. A check of the supposed dimensional reduction from four to 
three dimensions is of particular interest in connection with the electroweak phase 
transition. $ 6 terms in the dimensionally reduced S77(2)-Higgs model are usually 
argued to be irrelevant even at the transition to the spontaneously broken phase, 
and hence dropped pl| . In our extension of the previous investigations we will keep 



the $ 6 term in a 4d effective scalar theory at finite temperature, which is derived 
from the underlying S£/(2)-Higgs model by integrating out the gauge field degrees 
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of freedom. The phase structure of the effective scalar theory will then be studied 
in a finite and infinite volume. Hopping parameter expansions are supposed to work 
the better the smaller k, thus the larger the Higgs masses. Hence this investigation 
complements the range of Higgs masses which has been available in recent Monte 
Carlo simulations ||22|| . One of our aims is to find the critical Higgs mass above 



which the electroweak phase transition ceases to exist. (In case that the physical 
Higgs mass lies above the critical Higgs mass, it is bad news for an explanation of the 
observed baryon number asymmetry in the universe. The necessary ingredient for an 
out-of-equilibrium situation can no longer be provided by the electroweak transition, 
if the "transition" turns out to be truly a smooth crossover phenomenon.) 

A further application of our computational tools are (tri) critical phenomena in 
statistical physics. The order of the transition in a superconductor of type II has 
been recently under debate ( ||23|| and references therein.) The existence of a tricritical 
point for a suitable Ginzburg-Landau parameter has been conjectured [B3J], but a 
proof of its existence is still outstanding. Work in both directions is in progress. 

A Large coupling limit of O(N) lattice models 

For JV > 1, on a D-dimensional hypercubic lattice A, we consider the partition 
function 



Z = J £>$exp(-S($)) 
with corresponding expectation values 

P >x,a= \ J D$P($) exp (-£($)), (47) 



< 



where 



and 



2 

x^y 



S($) = $^ + A($ 2 -l) 2 + a($ 2 -P 3 

Measure and action are globally O(N) invariant. The observable P should be ap- 
propriately bounded so that the integrals exist. Fields at different lattice sites x 
and y interact by the hopping coupling v xy , which is assumed to obey the following 
conditions 

v X y = v yx = v(x-y), 
v(0) = 0, 
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^Vxy < OO, (48) 

Yl ''n^'r 'I',) 2 > 0. 

As an example, these conditions are satisfied for the pure nearest neighbour inter- 
action, where 

D-l 

V x y = 2K} j (S Xt y+p, + 6 X ,y--j{), (49) 

with ju denoting the unit vector in /^-direction. We consider fl47|) in the limit a —>■ oo 
with A = aa and a a fixed real number, i.e. 

< P > a = lim < P > aa>a 

for 

S ($) = $ 2 + a($ 2 - 1) 2 ($ 2 - 1 + a). (50) 

The hopping parameters v xy are considered as fixed. We have to select all field 
configurations that minimize the action in this limit. It is convenient to parametrize 
the fields according to 

§x = u x v x , u x e SW-i (the N-l sphere), v x > 0. 



Lemma A.l As a —>■ oo we get the following behaviour of expectation values P in 
dependence on a 

a > 1 : 

< P >a= K II f / d^N-i{u x ) J P(u) exp (-S(u)), 

S(u) = —- ^ v xyU x ■ U y . 
x,y£A 

a = 1 : 

< p > 1= Mi JJ I ( / dMN-l(Ux)Sv*,l + K,o) ) p ( uv ) ex P (-<%, 

xEA W=0,1 / 

5(tt, V) = - i ^ V X yV X VyU X ■ Uy. 

xgA x,ydA 

-3 < Q < 1 .- 

<P> Q =P(0). 
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a = —3: 



x)fiv x ,l + fiv x ,o) 

x€A \v x =0,l J s N~i 

P(V3uv) exp (—S(u, v)), 

S(U,V) = ^3vl - - ^ v xyV x V y U x ■ Uy. 
ieGA x,y£A 



a < — 3: 




< P >a= Na | | ( / dfitf-lK) ) P((l " ^)U) exp 



3 

2 ( x ~ Y^ VxyUx ' My ' 

x,j/GA 

The J\f a are positive normalization factors independent of P such that < 1 > Q = 1, 
<if2jv_;L £/ie standard measure on the sphere. 

For a > 1 and a < —3 we obtain the O(N) Heisenberg model (Ising model for 
N=l). For — 3 < a < 1 the lattice model becomes completely decoupled. At the 
boundary points a = 1 and a = —3 the result are "diluted" O(N) models, i.e. O(N) 
models with additional occupation number variables v x G {0, 1}. 

Outline of the proof: The properties (48) ensure that the minimizing field con- 
figurations of the action 5($) are translation invariant. Thus it is sufficient to 
determine the minima of 

F($) = ($ 2 - 1) 2 ($ 2 - 1 + a), $ G H N 

for the various values of a. Finally the saddle point expansion for the different cases 
yields the lemma. 

An alternative way to study the large coupling limit is to perform it termwise 
in the HPE series of correlation functions. We specialize to the nearest neighbour 
interaction (fl"9f). The only way the coupling constants A and a enter the linked 

cluster expansion is via the connected one-point vertex couplings v 2n (A, a), defined 
in Eqn. (|TTJ) . The connected one-point vertex couplings are related to the full one- 

O 

point couplings v 2n , defined by 

fd N ®$l n exp (-<?($)) , K1 . 

V 2n = 5 , (51) 

/ d N $exp(- S ($)) 

by the identity ([12]). At any finite order /, the coefficient of k is a polynomial in 

O c O 

the v 2nJ hence in the v 2n - Invoking a saddle point expansion again, we obtain 
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Lemma A. 2 Let X = aa with a > 0, — oo < a < oo. Define for nonnegative 
integers N, k 

^- (2fc - 1)!! r#fe' 

As u — ► oo we get for every k > 
a > I : 

V2k = A N , k + Oia~h. 



a 



-3 < a < 1: 



a = — 3: 



. 4- + 0(«7-s), , forN = l, 
l + 0(a"2) , forN>2. 



(2fc-l)!! 1 „ ~, 

v 2fc = i —L (i + 0(a 2)). 



3* 



T + 0(«7-5), , forN 



a < — 3: 



v 2k = A N k- < +1 , 

3 fc + 0((x-5) , forN>2. 



^2fc = A Njk (1 - y)^ + 0{<T-$\ 



For a > 1 the vertices are identical with those of the O(N) Heisenberg models. 
In the range — 3 < a < 1 they agree with the vertices of a purely Gaussian model 
with ultralocal action 

S ($) = (3-2a)<r$ 2 , 

leading to complete disorder as o — > oo for every finite k. For a < —3 we obtain 
the Heisenberg model again, as can be seen as follows. For any (5 > 0, rescaling of 
the vertices 

v 2 *- P 2k v 2k (52) 

implies a corresponding rescaling of the connected vertices, v 2k — > (3 2k v 2k , cf. (0). 
In turn, elementary graph theory shows that all susceptibilities change according to 
Xn(ft) — > fi n Xn{fi 2 K)- Hence, universality classes are invariant under (|52|) . Further- 
more we see that, for N > 2, the boundary points a = 1 and a = —3 belong to the 
Heisenberg class as well. A remnant of the occupation number variables is only seen 
in the case of N — 1, which is a remarkable exception and needs further study. 

Thermodynamic quantities like \ and the critical coupling k c can be directly 
determined in these limiting models. Alternatively, one may start with the original 
action (0) at finite a, calculate x an d ^ c in the HPE, and take the large coupling 
limit last. Our results indicate that both limits commute. 
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